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ABSTRACT 

In order to find more examples of the elusive high-redshift molecular absorbers, 
we have embarked on a systematic discovery program for highly obscured, radio-loud 
"invisible AGN" using the VLA Faint Images of the Radio Sky at Twenty centime- 
ters (FIRST) radio survey in conjunction with Sloan Digital Sky Survey (SDSS) to 
identify 82 strong (> 300 mJy) radio sources positionally coincident with late- type, 
presumably gas-rich galaxies. In this first paper, the basic properties of this sample 
are described including the selection process and the analysis of the spectral-energy- 
distributions (SEDs) derived from the optical (SDSS) + near-IR (NIR) photometry 
obtained by us at the Apache Point Observatory 3.5m. The NIR images confirm the 
late-type galaxy morphologies found by SDSS for these sources in all but a few (6 of 
70) cases (12 previously well-studied or misclassified sources were culled). Among 70 
sources in the final sample, 33 show galaxy type SEDs, 17 have galaxy components to 
their SEDs, and 20 have quasar power-law continua. At least 9 sources with galaxy 
SEDs have -fT-band flux densities too faint to be giant ellipticals if placed at their pho- 
tometric redshifts. Photometric redshifts for this sample are analyzed and found to 
be too inaccurate for an efficient radio-frequency absorption line search; spectroscopic 
redshifts are required. A few new spectroscopic redshifts for these sources are presented 
here but more will be needed to make significant progress in this field. Subsequent 
papers will describe the radio continuum properties of the sample and the search for 
redshifted H I 21 cm absorption. 

Subject headings: galaxies: photometry — infrared: galaxies — quasars: absorption 
lines 
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1. Introduction 

While most strong radio sources are identified with luminous elliptical galaxies or optical/UV 
point sources (quasars), a small number of these sources lie behind obscuring screens of gas and 
dust becoming invisible to us at optical/UV and possibly also at near-IR wavelengths. Identifying 
highly obscured, "invisible AGN" can lead to the detection and study of highly-redshifted radio- 
frequency atomic and molecular absorption line systems, which can provide important information 
about the early universe unobtainable in other ways. Radio absorption lines provide unparalleled 
velocity resolution and detailed physical diagnostics. And, where the background radio source is 
spatially-resolved, very- long baseline interferometry maps in the H I absorption line can provide 
the best spatial resolution (~ 10 pc at 2=1-2) of the gas distribution in distant galaxies; magnetic 
field st ructures also ca n be mapped using Faraday rotation measures of a polarized background 



source (|Gaenslerl 120071 ) . or using Zeeman splitting in the 21 cm line itself in the future. Atomic 
plus molecular absorption can provide much additional information. First, molecules trace the gas 
that supplies star formation so that the strengths and ratios of their absorption lines can be used 
to determine the physical conditions (e.g., temperature, pressure, density, chemical abundance and 
magnetic field strength) of the locations where most stars have formed over cosmic time. Second, 
the combination of various H I and molecular absorption/emission lines makes it possible to probe 
whether and by how much dimensionless fundamental physical constants have changed as the 
universe evolves. 

So far there are about 80 H I 21 cm absorbers detected at z > 0.1. Only five of them are 
molecular absorber^. The first system, PKS 1413+135, was detected in H I absorption in 1992 
at z = 0.25 (Carilli, Permian & Stocke 1992), in CO absorption in 1994 (Wiklind & Combes) and 
subsequently in nearly a dozen species (Wiklind & Combes 1997). This source is a BL Lac object 
residing in an edge-on spiral galaxy and the nucleus is highly obscured (Ay > 20; Stocke et al.1992; 
Carilli et al. 1992). Then B3 1504+377 was found at z = 0.67 (Wiklind & Combes 1996a), which 
is also an intrinsic molecular absorption system. The other three molecular absorbers are all in 
gravitational lens systems: B0218+357 at z cra = 0.94 with the lensing object at z^ s = 0.69 (Wiklind 
& Combes 1995, Menten & Reid 1996, Gerin et al. 1997), PKS 1830-211 at z em = 2.5 with the 
lensing object at z abs = 0.89 (Wiklind & Combes 1996b) and PMN J0134-0931 at z cm = 2.2 with the 
lensing object at z a b s = 0.77 (Kanekar et al. 2005, 2012). OH lines are detected in a ll five systems 



but only two have co njugate OH molecu lar absorption /emissi on lines (PKS 1413+135. lDarling||2004 



work ( Darlind 200 



Kanekar et al. 2004 : PMN J0134-0931. iKanekar et al. 1120121 ) important for fundamental constants 



Although "blind" searches of samples of strong radio sources have been undertaken to find 



1 In this paper molecular absorption refers to the multitude of species seen in the radio-frequency spectra of giant 
molecular clouds, not the UV absorption transitions of H2 and CO absorption arising in cool halo gas probed by 
far-UV spectroscopy in the Milky Way and high- z QSO spectra of intervening galaxies ( Noterdaeme et alj|20ld . 



Srianand et aPbood ) 
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molecular absorption, mostl y they have led to negative results, alth ough a few new H I 21 cm 
absorbers have been found (jCurran et al.l l2011al . iGupta et al.l 120091 ). Some of the unsuccessful 
surveys for high-z molecular absorbers include: millimeter observations of strong sources with 
known redshifts (Willett, Kanekar, & Carilli, private communication), cm- and m-wave observations 
of reddened quasars and radio g alaxies, sources with dam ped Lya absorbers, and type-2 quasars, 
etc. (Curran et al. 2006, 2011; ICurran Sz Whiting! l2010l ). While a few new H I 21 cm absorbers 
have been discovered, no new molecular absorbers have been discovered in these searches. This 
suggests both that high-z molecular absorbers are very rare and that rather standard techniques 
cannot be used to discover molecu lar absorbers (e.g., searches using samples containing damped 
Ly a absorbers; IGupta et al.ll2009l ); i.e., the obscuration we are trying to find significantly absorbs 
and reddens the quasar optical- UV continuum so that it is not detectable in most cases. Thus, 
optically-br ight quasars are not a go od sample to use to search for radio-frequency molecular 



absorption (jCurran et al. 



2008 



2011b! ). 



With the overall goal of discovering new examples of high-z molecular absorbers, we have 
embarked on an extensive discovery program using a new selection method devised specifically 
for this task. In this paper we describe the first steps towards the ultimate goal of discovering 
a substantial sample of high-z radio absorption line systems. This paper is arranged as follows: 
In Section 2 we describe this new selection method which co mbines Very Large Array (VLA) 
Faint Images of the Radio Sky at Twenty centi meters (FIRST: iBecker et al.l Il995l ) radio source 
detections with Sloan Digital Sky Survey (SPSS: lYork et al.ll2000l ) photometry and morphologies 
for the coincident galaxies designed to select strong radio sources which are associated with late- 
type galaxies. In Section 3 we describe new broad-band near-IR (NIR) imaging of this sample 
obtained by us at the Apache Point Observatory's 3.5m telescope (APO hereafter) which allows us 
to search for nuclear point sources which are obscured optically and to extend the spectral energy 
distributions (SEDs) of these sources from the SDSS optical photometry to longer wavelengths. 
The analysis of the NIR images and of the optical/NIR SEDs is described in Section 4. A few 
new spectroscopic redshifts obtained in the progress of this work are also reported in this Section. 
Individual results for some representative and also unusual objects are presented in Section 5 and 
candidates for finding atomic and molecular absorption are identified. We conclude with a summary 
of the observational results and discuss the future work using our on-going survey in Section 6. 

Paper 2 in this series will describe high-resolution VLA and Very Long Baseline Array (VLBA) 
radio maps of these sources to identify those sources containing very compact structures and to 
identify sources which have radio/optical-NIR position offsets indicative of foreground, not asso- 
ciated absorption. A first round of H I 21 cm absorption observations from this sample will be 
presented in Paper 3. 
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2. Sample Selection 



2.1. Background 

The host galaxies of radio-loud AGN are almost exclusively giant ellipticals (e.g. Urry &; 
Padovani 1995, McLeod 2006). This is true for all the subclasses of radio- louds including quasars, 
FRI and II type radio galaxies and BL Lac Objects. While we might generically expect radio sources 
associated with a spiral or irregular host galaxy or with a merger to possess foreground absorption, 
their occurrence is extremely rare. But recent large surveys allow us to find rare objects that do 
not obey the general tendency for radio-loud AGN to be hosted by giant ellipticals for specific 
reasons. For example, Wilson &: Colbert (1995) have suggested that the origin of radio-loud AGN 
is in the merger of two supermassive Black Holes (SMBHs) when two disk systems merge to form 
an elliptical galaxy (Mihos &: Hernquist 1996). If the Wilson & Colbert (1995) scenario is correct, 
the merger of the SMBHs first creates a luminous but small-scale (~ lkpc) r adio-loud AGN termed 



a "Compact Symmetric Object" (CSO; lBegelmanlll996l ; iPerlman et al.ll200ll ) which can then evolve 
into a large-scale classical double (FR II) through subsequent outbursts. And if the CSO forms 
quite quickly after t he merger, it can be accompanied by remnant gas and dust from the merger 



(jWillett et al.ll2010l ). This circumstance can cause the galaxy's optical morphology to appear 



distorted and provide gas for fo reground absorption at the same time. Indeed, some CSOs have 
disturbed optical morphologies (IPerlman et al,ll200ll ) and a high fraction (^ 20-40%) of compact 



radio sources including CSOs have H I absorption in their radi o spectra (Ivan Gorkom et a. 



Vermeulen et al.l 120031 : iGupta et al.l 120061 : IChandola et al.l l201ll : but see ICurran Whiting 



19891 : 
2010). 



PKS 1413+135, the first molecular absorber detected, is also a CSO ( Perlman et al.lll996l ) with 
intrinsic absorption. 

Besides the two intrinsic molecular absorbers, known gravitational lenses have provided 3 other 
known molecular absorbers, B2 0218+357, PKS 1830-211 and PMN J0134-0934. Gravitational 
lenses possess unusual optical morphology due both to their multiple images and to the presence 
of a foreground galaxy. It is also possible that a foreground, gas-rich galaxy can produce atomic 
and/or molecular absorption even without producing multiple images (Narayan & Schneider 1990). 
Again, a correlation between unusual optical morphology and foreground absorption can result, 
although in this case the absorbing gas is external to the AGN and its host. 

Another rare class of radio-loud AGN which possesses unusual optical morphology is the so- 
called "alignment effect" radio galaxy (2 known examples are in our sample) in which UV emission 
lines from the gas surrounding the radio jets are redshifted to the optical and are responsible for a 
non-elliptical morphology in certain wavebands (e.g., Zirm, Dickinson & Dey 2003). This circum- 
stance has not led to the discovery of molecular absorbers, probably due to the very extended radio 
source structure required to create the "alignment effect". On the other hand, the first two classes 
described above often include radio sources with compact radio structure. Compact structure fa- 
cilitates the detection of atomic and molecular absorption because the absorption detection does 
not get diluted by flux which does not go through the absorbing gas (see Paper 2). 
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2.2. The Radio-loud Late-type Galaxy Sample 

The 1.4 GHz VLA FIRST survey covers over 9,000 square degrees of the North and South 
Galactic Caps. FIRST has a typical flux density rms of 0.15 mJy, and a resolution of 5". The Data 
Release 5 (DR5) of SDSS covers about 8,000 square degrees of almost exactly the same area of 
the sky as FIRST. This field includes approximately 1500 FIRST sources with flux densities > 300 
mJy. Among these sources, 808 are associated with SDSS sources, i.e., one or more optical objects 
are detected within 1.5" of the radio centroid; thus, it is possible that some of the associated SDSS 
sources are not the AGN host galaxy. Excluding 467 sources classified as stars which might actually 
be QSOs, 341 are probably the host galaxies of the radio sources . Note that star-galaxy separation 



is robust down to r ~ 21.5 mag in the SDSS (jLupton et al.ll200ll ). To select non-elliptical galaxies 



two standards are used. The first is that the best-fit of the radial profile is exponential (typical for 
spiral galaxies) rather than De Vaucouleurs (typical for elliptical galaxies; 94 sources). The other 
is that the core compactness parameteB is smaller than 2.6 (also typical for spiral galax.es; 200 
sources). These two criteria also eliminate quasars and point-source (AGN) dominated galaxies 
as well. The combination of the two standards makes the selection fairly reliable (Strateva et al. 
2001) but there could still be some mis-classified cases. All SDSS classifications employed the i- 
band images (rest-frame g-band at the median photometric redshift of the sample, z ~ 0.6), which 
typically have the highest signal-to-noise images for our sample. 

Overall, 82 objects in SDSS DR5 meet both these criteria (see Table 1). None of these sources 
has more than one optical counterpart, although the possibility remains that some of the optical 
sources are actually merged images of two or more galaxies. Because we were interested in ob- 
taining a representative sample of these rare AGN, we did not conduct any completeness tests for 
this sample, nor did we obtain a new sample using later SDSS data releases. Two nearby galax- 
ies are excluded because their intrinsic radio luminosities are much lesser than typical radio-loud 
AGN (J1413-0312 = NGC 5506 and J1352+3126 = UGC 8782). Also we removed J0834+1700 
because the SDSS galaxy is coincident with one of the two double lobes in a classical double 
radio source; an unobscured elliptical galaxy is located between the double lobes, and so is the 
more likely counterpart to the radio source. We have also excluded 4 already well-studied ob- 
jects: PKS 1413+135 (J1415+1320) mentioned above, another well-known low-z CSO 4C +12.50 
(J1347+1217), and the "alignment effect" radio galaxies, 4C +23.27 (J1120+2327, @ z=1.819) 
and 4C+24.28 (J1348+2415, @ z=2.879), which lack core radio emission. Our sample selection is 
demonstrably successful because all the three possible situations mentioned above have a few known 
examples in the sample. PKS 1413+135 is the first discovered high-z molecular absorber and a 
CSO. B2 0748+27 (J0751+2716) is a gravitational lens (Lehar et al. 1997). And 4C +23.27 is the 
well-known "alignment effect" radio galaxy 3C 256 (Zirm, Dickinson & Dey, 2003) just mentioned. 
There are a couple of other known CSOs and Compact Steep Spectrum Sources (CSSs; slightly 



2 The r atio of the Petrosian r adii containing 90% and 50% of the galaxy light. Same as the the concentration index 
defined in lStrateva et al.1 ([200 ll ). 
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larger in physical extent of 1-20 kpc but similar to CSOs; Ide Vries et al.lll998l ) present and more 
have been found in the course of this work (Paper 2). A few host galaxies appear to be merging 
with another galaxy. 

The full list of all 82 sources in our sample is found in Table 1. This Table includes by columns: 
(1) coordinate style designation for the source which is used in this paper; (2) the commonly used 
radio source name for the source; (3) the 20 cm continuum flux density from FIRST in Jy; (4) the 
associated SDSS name for the optical counterpart which also provides an accurat e (RA,DEC) cen 



troid position; (5) the SDSS r-band "model magnitude" and its associated error (jStoughton et al 



2002J) , which results from an extrapolation of the best-fit (for this sample) exponential radial pro- 
file; (6) the photometric redshift and error from Data Release 8 (DR8) using an empirical method 
called "kd-tree nearest neighbor fit" (Csabai et al. 2007); and (7) the spectroscopic redshift where 
available. The spectroscopic redshifts come from references in the NASA/NED database except 
for the cases noted with an asterisk; these new spectroscopic redshifts were obtained by us and are 
described in Section 3. 

Notice that most of these sources are quite faint in the optical (~ 75% fainter than r=21) so 
that the application of the Strateva et al. (2001) morphology discriminators is not without some 
uncertainty. For this and other reasons we sought to confirm the diffuse nature of the potential 
invisible AGN hosts using more recent SDSS results (DR8) and new near-IR imaging (see Section 
3). Since the original galaxy classification was made using SDSS DR5, we have revisited the sample 
selection process using DR8. In almost all cases the original classification is confirmed but in a 
few cases (5 in number) the DR8 compactness parameter (C) exceeds the C < 2.60 cut suggested 
by Strateva et al. (2001) for late-type systems. In one other case (J1421— 0246) the r-band SDSS 
galaxy image is much better fit by a De Vaucouleurs r 1//4 -law than an exponential disk. Our near-IR 
imaging of these sources confirms these reclassifications, except in one case. For J0824+5413, our 
near-IR imaging reveals a very diffuse galaxy, at odds with the C=2.9 compactness value for this 
object in DR8. Therefore, we have retained this source in our sample. The other reclassifications 
are marked with daggers (f) in Table 1 and removed from the final sample listed in Table 2. 



3. Observations 

3.1. Near-IR Imaging 

In order to check the SDSS galaxy classifications, to complement the SDSS photometry (and 
thus to provide an optical/NIR SED), and to search for an extincted point-source indicative of an 
imbedded AGN, we obtained NIR photometry of this sample using the Near-Infrared Camera & 
Fabry-Perot Spectrometer (NIC-FPS; Hearty et al. 2005) in its imaging mode on the APO 3.5 m. 
In the last several years almost all of the objects in this sample were observed in the J and -fCjhort 
bands with all objects detected in -fT s hort and most in J-band. ii-band observations were made for 
fewer of these sources, particularly where the ff-band photometry was important in evaluating the 
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nature of the source. During these observations, the seeing varied from sub-l" to 2" with a median 
seeing value of ~ 1.3". This image quality is approximately the same or, in many cases, better than 
the image quality for the SDSS. 

The images were reduced in the standard manner with flat-field and dark images obtained 
the same night. Because all of these images were obtained using a 5-point dither pattern on the 
sky, the science images were used to create a sky-flat by median-filtering. Following flat-fielding 
the images were aligned and stacked; astrometric positions for source centroids were set to the 
World Coordinate System using 5-10 SDSS stars identified in these images. Most of the objects 
in this project are faint, so, in order to reduce the photometric error and for consistency with 
the SDSS optical photometry, we first measure the magnitude in an aperture of 2" (my)- This 
magnitude is similar to the fiber magnitude maber in SDSS (Stoughton et al. 2002). SDSS also gives 
the "model magnitude" m mo dei as the best estimated magnitude for extended sources. Assuming 
an unchanged radial profile between the optical and NIR images, the "model NIR magnitude" is 
m = my + (w-modei — ™fiber) and is the NIR equivalent of the SDSS model magnitudes to facilitate 
a consistent source SED. 

These are the magnitudes that we quote in Table 2 and throughout this paper. Table 2 lists 
our NIC-FPS photometry for the 70 sources in the trimmed sample including by column: (1) 
source name as in Table 1; (2-4) the next three columns list the J, H and -Kshort band magnitudes. 
The listed statistical errors include a contribution from photon statistics as well as the zero-point 
errors from the 5-10 2MASS stars in the field. Column (5) lists the UT date of the observations in 
the format YY-MM-DD; all observations of any one source were obtained sequentially on the UT 
date listed. Column (6) lists the classification of the source based on its optical/NIR SED. These 
classifications are discussed in the next Section in detail. Abbreviations include: G = galaxy; Q = 
quasar and Q + abs = possible reddened quasar; G + Q = galaxy + quasar. Column (7) lists zg t , 
the best-fit photometric redshifts based on optical + NIR SEDs for G-type sources only. Column 
(8) lists the Hubble type for the best-fitting SED galaxy template (see next Section for discussion) 
used for the classification in Column (7). The entries in the last three columns will be discussed in 
detail in the next Section. 

Usually there are 5-10 non-saturated objects from the Two Micron All-Sky Survey (2MASS) 
catalog in the ~ 5' x 5' NIC-FPS field, which are used to provide zero-points for the magnitudes 
listed in Table 2. The magnitude error budget includes amounts from the image photon-statistics 
(usually o"2» < 0.05), the 2MASS magnitude system error (a < 0.1) and slightly larger and less 
well-defined errors accrued in placing the NIR magnitudes onto the SDSS "model magnitude" scale 
(cmodci ~ 0.1 and o"fib er varies from 0.05 to 0.3). Sources which show NIR point sources, and so 
have more centrally-concentrated surface brightness profiles than their optical images, may have 
NIR "model magnitudes" which have been under-estimated systematically by 0.1-0.2 mags. For 
non-detected objects, we set a 3a flux level recorded as a model magnitude limit in Table 2. Blank 
spaces in Table 2 mean that the source was not observed in that band. 
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The systematic astrometric error of an SDSS frame is on the order of 50 mas per coordinate. 
For objects brighter than r = 20 mag, the object controiding for a point source gives a random 
error of about 20 mas. About 5-10 SDSS stars with r < 20 mag were selected in each NIR field for 
astrometric calibration, yielding a total positional uncertainty in the range of 50 - 100 mas. For 
bright FIRST radio sources the astrometric error is smaller than 50 mas so that the optical/NIR 
positional errors dominate the accuracy with which the SDSS, NIR and VLA/VLBA images can 
be registered. All sources have coincident NIR and optical positions although discrepancy exists 
between optical/NIR and radio positions for some of the sources mostly due to extended structure 
of the radio sources. More details will be discussed in Paper 2. 

Specific unusual circumstances involving morphology, relative positional offsets and unusual or 
inconclusive SEDs are discussed source by source in Section 5. 

Due to the faint and diffuse nature of most of the target galaxies in this sample, spectroscopic 
redshifts are available for only 17 of them. Originally it was thought that a reliance on the SDSS 
photometric redshifts would be sufficient for this program given the wide bandwidths available on 
current generation sub-GHz radio receivers. However, this has not proven to be the case. First, the 
special defining qualities of this sample make photometric redshifts (photo- z hereafter) less reliable 
than for most galaxy samples (i.e., optical objects can be combinations of starlight plus extincted 
AGN). And second, the difficult radio frequency interference (RFI) environments now present at 
all terrestrial radio telescope sites make even typical photo- z accuracy insufficient for a sensitive 
H I and molecular line search; i.e., without an accurate spectroscopic redshift, the presence of 
intense, receiver-saturating RFI means that it is not possible to determine definitively whether an 
absorption is present without a specific frequency prediction. Due to this circumstance we initiated 
a program for obtaining optical and/or NIR redshifts using the spectrographs at the APO 3.5m. 



3.2. Optical and Near-IR Spectroscopy 



The Dual Imaging Spectrograph (DIS II) was used in its standard low dispersion mode to obtain 
150 km s" 1 resolution spectra in the 4000-9000 A wavele ngth region (1.8 A pix shortward of 



5500 A and 2.4 A pix longward; see lOke Gunnl Il982l for original design). TripleSpec is a 
cross-dispersed NIR spectrograph that provides simultaneous R ~ 3500 spectra in a continuous 
wavelength co verage from . 95-2. 46 \xm in five spectral orders. The instrument is described in 
more detail in lwilson et al.1 ( 2004 ). Where we have successfully obtained DIS II spectra, emission 
lines are detected and redshift errors are ±30-60 km s _1 and ±90 km s _1 where absorption and only 
weak emission is detected. TripleSpec NIR spectra have similar resolution but the signal-to-noise 
ratio is much lower so that the DIS II redshifts are preferred where we have obtained both optical 
and NIR spectroscopy. These new redshifts led to the detection of redshifted H I 21 cm in two 
cases. Despite observing many more objects which possessed only photometric redshifts only one 
new redshifted H I 21 cm absorption detection was made based on a photo-z alone: J1357±0046 
(PKS 1355+01) possesses two H I components at zhi=(0.7971, 0.7962); i.e., at a significantly higher 



- 9 - 



z than the SDSS z p h o t = 0-57. The eight new redshifts we have obtained for this sample are at the 
limit of what is possible using a 4m class telescope; i.e., r < 20 and/or K < 16.5. The spectra 
we have obtained thus far for this project using DIS II and TripleSpec are listed in Table 3 with 
the instrument used, the redshift obtained including measurement errors and the emission lines 
detected. These spectra are shown in Figure [Q and Figure [2j 

4. NIR Morphology, Magnitudes and Optical/NIR SEDs 

With multi-band observations, we are now able to better understand these sources despite 
their faintness. In this Section we use the NIR magnitudes to examine the NIR morphology of the 
source, use the NIR magnitudes to create an optical-NIR SED and compare these with template 
galaxy spectra. The -KTshort magnitudes are then used to place these objects on the K — z diagram 
for which many radio-loud sources show a strong correlation. 

4.1. Near-IR (NIR) Morphology 

The NIC-FPS ( JHK) imaging provides an independent check of the SDSS galaxy morphology 
classifier. In general, these images provide ample confirmation of the galaxy morphology in that 
virtually all of them are diffuse sources in moderate seeing conditions (median ~ 1.3", Figure [3]). 
For a few sources, particularly some of the optically red sources, the NIR images suggest that an 
elliptical galaxy morphology is possible so that the Strateva et al. (2001) classification may be 
incorrect. The few sources which show a centrally-concentrated structure in our NIR images are 
J1354+5650, J0749+2129, J0920+2714, J1043+0537, J1019+4408,and J0003-1053. These are not 
removed from the sample because they may still be extincted, although this is less likely compared 
to the more diffuse galaxies (e.g., J0920+2714 is detected in redshifted H I 21 cm absorption). 
Otherwise, the optical morphology classifications for these sources made by the SDSS is confirmed. 

4.2. Optical-NIR Spectral Energy Distributions (SEDs) 

These sources have a wide variety of optical-NIR SEDs due to three factors: variations in host 
galaxy type, the relative importance of the host galaxy and the AGN as a function of wavelength, 
and the amount of nuclear obscuration. In an attempt to categorize their SEDs we have placed each 
source into one of the following three categories, listed these in Table 2 and show a few examples 
of each type in Figures HHS]) : 

• Quasar (Q in column 6 of Table 2) or obscured quasar (Q+abs). The Q SED is approximated 
by a single power-law declining from blue to red in Fa, although in some cases contributions 
from emission lines to specific flux bands are evident. Almost all previously known quasars in 
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our sample have Q-type SEDs. Sources classified as Q+abs show evidence for a flattening in 
flux at the blue end. In some cases the SDSS u-band flux is only an upper limit making the 
classification uncertain between Q-\-abs and a galaxy G-type SED; i.e., it is uncertain whether 
the near-UV flux continues to decline or not. See Figured] for examples. 

• Galaxy (G). This galaxy-like SED peaks in the middle bands and declines toward both the 
blue and the red ends. See Figure [5] for examples. If the turn-down at the blue-end is not 
explicitly observed (u and sometimes g-band fluxes are upper limits), the source SED could 
either be Q+abs or G. Sources are classified as Q+abs if the NIR slope is too steep to be a 
galaxy, otherwise these sources are classified as G. 

• Galaxy + Quasar (G+Q). The SED is a power-law Q component with a bump in the middle 
bands which could be a galaxy G component. Sources missing the flattening at the red end 
are labeled as G+Q since the absence of a flattening or turn-down in the NIR could be due 
to contributions from a nuclear source that is obscured at shorter wavelengths. See Figure [6] 
for an example. 

Our classifications are largely confirmed by the 17 sources for which optical and/or near-IR 
spectroscopy is in-hand. Overall there are 33 sources in the G class (including 12 SEDs best-fit by 
early-type galaxies and 21 by late-type galaxies), 17 in G + Q, 6 in Q and 14 in Q + abs. While 
our NIR photometry was obtained several years after the SDSS optical photometry, we do not see 
strong effects in the SEDs arising from any AGN variability that may have occurred over that time 
period. The SDSS and NIR images (as in Figure 3) and the SEDs (as in Figure 4-6) will be shown 
for the entire sample in Yan (2012, PhD thesis, in preparation). Due to the strong possibility that 
the 33 sources in the G class are actually rather simple in ther overall SEDs, in this next Section 
we concentrate on G-class sources in a further scrutiny of their photometric properties. 



4.3. Photometric and Spectroscopic Redshifts 



Since up to now this program has relied primarily on the SDSS photo-zs, a brief review of their 
basic properties is in order. Before Data Release 7 (DR7), S PSS photo-zs are der i ved from fits 



using typical SDSS resolution spectra of various galaxy types (IBudavari et alJl2000l ; ICsabai et al 



20031 ). A new algorithm is used in DR7 which utilizes SDSS spectros copic catalog as a training set 
to determine photo-zs from spectroscopic objects with similar SEDs (jCsabai et al.ll2007l ). However, 
objects in this training set are brighter and closer on average than typical objects in our sample. 
A fainter and deeper training set more representative of our sources is included in DR8 which 
has r-band magnitudes mostly between r = 20-22, redshifts mostly between 0.2-0.8, and the 
improved average error is 0.10-0.12 (Csabai, private communication). Although these redshifts are 
a good potential match to our own sample, the correlation between the spectroscopic and photo- zs 
we have in-hand is less compelling. Figure [7^ shows sources with known spectroscopic redshifts 
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compared with their photometric redshifts from DR8. Only half (8 out of 16) of these sources 
appear reasonably well-fit by the DR8 algorithm to the la accuracy noted above. But even for the 
8 with Zphot — z spec , the photometric redshift systematically underestimates the redshift by Az ~ 
0.05-0.1 (see Figure [7^,), i.e., the offset of H I absorption can be as large as 70 MHz. We note that 
the only source (J1414+4554) to the left of the solid line in Figure [Th, has an uncertain redshift 
of z = 0.186 (Falco et al. 1998). Given its unique and unusual location in this plot we doubt the 
accuracy of this spectroscopic redshift (Falco et al. did not publish their spectrum) and do not use 
it further. 

In the poorly- fit 7 sources at z > 1, three are classified as Q SED types and 4 are classified as 
Q + abs. It is understandable in these cases that the SDSS photo- z classifier obtains an unreliable 
redshift by trying to fit a galaxy template to a reddened quasar SED. An emission line spectrum 
with weak continuum may also pose a problem for the photo-z algorithms. Thus, those source 
with z spcc >> Zphot are not representative of the remainder of the sample and photo-zs for Q and 
Q + abs type sources in particular should not be used. 

Given the potential complexities of the SEDs for these galaxies, it would be very helpful if we 
could use our NIR photometry to further constrain source redshifts. But given the results shown 
in Figure [7H we restrict this analysis to only those sources classified as G-types in Figure 03. Using 
the SED templates in Figure [8] we have used a least-squares approach to determine the best-fit 
Hubble types and redshifts (zfit) for the G sources based on their SDSS + NIR SEDs and listed 
them in Table 2 column 7 & 8. We searched parameter space in redshift, peak flux density and 
discrete Hubble types (Figure [8]) to find the best-fit values that give the least difference between the 
observed flux points and the template SEDs. As shown in Figure [TJd using data on G-type sources 
only, a photo-z can be estimated to ± 0.2 at high confidence (~ 90%). While this is true both 
for those sources with well-fit SDSS + NIR SEDs, and for those sources with SDSS photo-zs only, 
there is a slight systematic difference between zg t and z p hot- This slight systematic underestimate 
of Az ~ 0.05-0.1 is the same as in Figure [7^, suggesting that the SDSS+NIR SEDs are a slightly 
better indication of the true redshift of these sources. However, we caution that this conclusion 
is based on a training set of < 10 sources. Even with this improvement in the precision, photo-zs 
are inadequate for a successful redshifted H I 21cm absorption since the above quoted precision 
(Az ~ 0.1) can provide a search frequency accurate to no better than ~ 40 MHz at z ~ 1. Given 
current levels of RFI at large radio telescopes world-wide this is still inadequate for a conclusive 
H I search (i.e., RFI- free regions are typically much smaller than this at v < 1 GHz) so that 
spectroscopic redshifts are still required. 

But before we can accept that the photo-zs for all G-type sources are accurate to ±0.1 in Az, 
we need to scrutinize them further using previous results for luminous radio galaxies. The G-type 
sources can be plotted on the well-known K — z diagram that relates the if -magnitude and the 
redshift (van Breugel et al., 1998). Willott et al. (2003) observed 49 radio galaxies with 151-MHz 
flux densities greater than 0.5 Jy and gave a K — z relation fitted by a second-order polynomial 
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between if -magnitude and log 10 z, 

if = 17.37 + 4.531og 10 z - 0.31(log 10 z) 2 . (1) 

The dispersion in the correlation between these two quantities for radio-loud AGN is reported 
to be 0.58 mag up to z = 3. If a galaxy is formed at z > 5, undergoes an instantaneous starburst, 
then evolves passively and has a luminosity of 3L* (typical giant elliptical galaxy hosting a radio- 
loud AGN) at present, the if -magnitude evolution of such a galaxy is close to the observed K — z 
relation presented here. The observed if — z values from our sample (excluding all Q and Q + abs 
type sources) are plotted in Figure [9] together with the K — z relation from Equation [TJ 

While some sample sources conform quite well to the relationship (including all those with 
spectroscopic redshifts), there are a significant number which should have z > 1 based on their 
if— band magnitudes and instead have z p hot < 1- While we will not understand these discrepancies 
fully until spectroscopic redshifts are obtained, there are several possibilities including: (1) the host 
galaxies of some radio-loud AGN are not giant ellipticals; (2) the host galaxy, or at least the inner 
bulge of the host, is heavily extincted (Ay > 20) even at if -band (unlikely); (3) the SDSS galaxy is 
not the AGN host but is rather a lensing, or simply a foreground, galaxy close to the sightline; and 
(4) the photo-zs are systematically too low due to the SDSS source being a composite of quasar 
and galaxy. Regardless of the solution for each case, the source offset from the K — z relation can 
be used as a diagnostic that there is some unusual circumstance for that particular source and that 
the photo- z is unreliable. This appears to be the case for a large fraction of the sample. 

Given Equation (1) and that the photo-zs are < 1, almost all G-type sources in Table 2 with 
if >17 would be discrepant in Figured! Approximately 2/3 of the sources falls into this magnitude 
range and so have zk >> 2fitj where zk is the redshift calculated in Equation (1) from the if -band 
mags in Table 2. The 21 discrepant objects in our sample have diffuse morphology in the NIR 
and so are either extincted even at if -band, are under luminous hosts, or are foreground galaxies 
unrelated to the AGN. Spitzer Space Observatory mid-IR photometry of these objects can determine 
if they are still extincted at if -band. In any case these sources are good candidates for foreground 
absorption. In Table 4 we list those 9 sources that are more than 2a fainter than expected in 
if -band based on their photo-zs. Table 4 contains the following information by column: (1) source 
name; (2) zg t from Table 2 column (7); (3) Hubble type from Table 2 column (8); (4) if -band 
magnitude from Table 2 column (4); (5) zk; and (6) the amount of redshift discrepancy, Az= 

ZK ~ Z&t- 



5. Discussion of Selected Unusual Sources 

J0003— 1053 This is one of the reddest sources in our sample with (r — if s) = 5.5 ± 0.3. 
Although the flattening of the SED longward of 1 micron is comparable to what is expected for 
a z ~ 2 elliptical, the observed blue flux density is much too faint for this interpretation. So this 
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source must be heavily obscured and is an excellent candidate for absorption studies at z = 1.474. 

J0134+0003 . The optical spectroscopic redshift of z=0.879 is from Drinkwater et al. 



(1997). ICurran et al.l (|2006l ) observed this source for H I absorption but severe RFI was present 



around the expected frequency for H I at z=0.879. 

J0751+2716 This source is a well-known "Einstein Ring" gravitational lens system with 
ziens=0.35 and zqso=3.2. While this source has been searched for redshifted H I 21 cm absorption 
by us, RFI was present which rendered the search inconclusive. 

J0901+0304 Our DIS II spectrum yields z = 0.2872 based on strong narrow emission lines 
(Table 3; Figure [JJ and led to the discovery of two components of redshifted H I 21 cm absorption 
at z=(0.2886, 0.2869). See Paper 3. 

J0920+2714 While this source is apparently associated with a rather normal bright 
(X=14.3) galaxy at z spcc = 0.2064 (see Table 3), the FIRST radio source position is offset from the 
SDSS galaxy position by 1.7 arcsecs, an offset we have confirmed using new VLA and VLBA maps 
(see Paper 2). An H I 21 cm absorption has been discovered at z= 0.2067 (see Paper 3), consistent 
with the optical/NIR galaxy; i.e., the radio-loud AGN is background to the galaxy which creates 
the absorption. 

J0939+0304 There is a point source seen only in the K s image of this source but not at the 
center of the "host" galaxy (Figure E} . Regardless of whether the AGN is in an interacting system 
or whether the nuclear offset is due to a foreground galaxy, this source is an excellent candidate for 
absorption studies. 

J1357+0046 The optical and NIR images show a faint linear string of knots, one of which 
is coincident with the radio source (see Paper 2). Despite the absence of a spectroscopic redshift, 
we have detected H I 21 cm redshifted absorption at z=(0.7971, 0.7962). See Paper 3 for details. 

J1502+3753 A 5 GHz MERLIN image shows that the very blue optica l/NIR counterpart 



sits in the center of a very extended, double- lobed source (jBolton et al.ll2006l ). It is an unlikely 
candidate for radio absorption studies because the radio source structure is so extended, with only 
a weak compact, core component. 



6. Discussion and Conclusions 

With the goal of finding more high-z molecular absorbers, we have selected 70 strong radio 
sources with non-elliptical morphology in the optical and obtained NIR images at the APO 3.5m 
telescope to further characterize their basic nature. Selection from the FIRST survey at /20cm > 300 
mJy assures that we are targeting radio-loud AGN. An SDSS associated galaxy classification as 
non-elliptical gives us the greatest probability of finding plentiful gas, regardless of the reason for 
the late- type morphologies; e.g., the observed galaxy could either be a very young radio galaxy like 
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a CSO or a foreground galaxy like in a gravitational lens system. In this first paper on "Invisible 
AGN" we have presented a selection method for finding highly-obscured AGN and the optical-NIR 
data necessary both to confirm their morphology and judge the amount of obscuration present. 

We have obtained NIR images of the sample to supply us with morphological information to 
check the SDSS classification, to extend the SDSS optical SEDs into the near-IR and to search for 
heavily obscured nuclei. Our NIR observations are consistent morphologically and photometrically 
with the expectation that many of these sources are heavily obscured systems. Analysis of the NIR 
morphology and of the optical/NIR SEDs provides a consistent picture of a heavily obscured, radio- 
loud AGN in most cases (see detailed discussion in Section 4.2). Also we have made a comparison 
between the photomeric redshifts and redshifts estimated from the -RT-band magnitudes. Because 
radio-loud AGN exhibit a significant "Hubble diagram" type correlation out to z ~3, we can use the 
K-hand magnitudes as a crude redshift indicator. For nine sources these K-band derived redshifts 
are significantly greater than their photo-zs (see Table 4 and Figure 7), indicating either that the 
host galaxy is not a giant elliptical or that the SDSS detected galaxy is forground to the AGN. In 
either case, these are excellent candidates for finding foreground absorption. However, this analysis 
also makes it clear that it is not advisable to use photometric redshifts of Q or Q+abs- type sources 
to conduct an efficient H I 21 cm absorption-line search. Additionally, Figure 7b shows that even 
when just the G-type sources are used, there is a systematic offset (Az ~ 0.05 — 0.1) between the 
SDSS-derived photometric redshifts and the spectroscopic redshifts on top of the standard statistical 
errors present in the SDSS photo- z measurement (Az ~ 0.1-0.12). While it may be possible to 
conduct an H I 21 cm search of sources with G-type SEDs using photometric redshifts derived using 
both SDSS+NIR photometry (called z&t herein), the statistical errors present (see Figure 7b) are 
too large to make such a search successful in the presence of the typical RFI environment found at 
ground-based radio telescopes. Again a more accurate spectroscopic redshift is required. 

While the optical/NIR properties described herein are compelling that most of the sources 
within sample are obscured, 'invisible AGN', it is also important to continue to scrutinize the basic 
properties of the radio sources against which absorption might be observed. Regardless of the 
amount of obscuration between us and the nucleus of an AGN, its radio spectrum is unlikely to 
show strong absorption lines if the radio source contains mostly geometrically-extended flux. This 
is because we are viewing the large majority of the source flux along unobscured lines-of-sight that 
avoid going through the host galaxy (and its obscured nucleus) entirely. VLA A-configuration and 
VLBA maps are in-hand and will provide identification of the most compact radio sources in this 
sample. These high resolution maps also will be used to make an accurate comparison between 
the radio source structure and the optical/NIR galaxy positions to identify potential intervening 
absorption candidates. These maps and their analysis will be presented in Paper 2. Not only will 
VLA/VLBA observations discover new examples of CSOs and CSSs but they will also identify those 
obscured nuclei which are likely to show the strongest absorption due to possessing very compact 
radio continuum structures. Hopefully, some of the absorptions we discover will be molecular lines. 
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Table 1. Basic Source Data. 
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i y 


4- 
ZIZ 


n 

u, 


01 

Ul 


n 
u 


90 
Z J 


4- 
ZIZ 


n 

u, 


04 

U4: 




700034-^01 9 

JUaUOn^OUl^ 


AC 4-^0 98 

4:V_/ T'JU.iO 


n 

u 




7000340 


81 4-^01 93^ 


o 


99 


1 

1 iJ 


4- 

ZIZ 


n 


99 




u 


40 

4U 


4- 

ZIZ 




u, 


99 




7000^4-41 98 


R3 nori94-4i r 


n 
u 


48 
-to 


7000^99 


1 84-41 9830 


c 

o 


91 
Z 1 


83 

,oo 


4- 
ZIZ 


n 

u, 


00 


n 
u 


U 


4- 
ZIZ 


n 

u, 


1 (\ 

1U 




juyu ( +0413 


4L- +04.30 


u 


( U 


J090 ( 50 


, (4+041oo( 


o 
O 


99 
ZZ 


,5 ( 


1 

± 


U, 


i n 
19 


U 


by 


1 

± 


0, 


1 1 

11 




J0910+2419 


4C +24.19 





.83 


J091022 


.55+241919 


,5 


20, 


.92 


± 


0, 


06 





,74 


± 


0, 


,14 




J0915+1018 


TXS 0912+105 





35 


J091512 


,95+101827 


,4 


22 


,01 


± 


0, 


,12 





,31 


± 


0, 


,03 




J0917+4725 


B3 0914+476 





35 


J091727 


,16+472525, 


,6 


21 


,82 


± 


0, 


,12 





,81 


± 


0, 


,22 




J0920+1753 


4C +18.29 


1. 


.08 


J092011 


.12+175324 


.7 


22 


.43 


± 


0, 


,16 





,70 


± 


0, 


,11 




J0920+2714 


B2 0917+27B 





.46 


J092045 


.05+271406, 


.6 


18, 


.04 


± 


0, 


,01 





17 


± 


0, 


,01 


0.206* 


J0939+0304 


PKS 0937+033 





47 


J093945 


.16+030426 


.4 


20, 


.44 


± 


0, 


05 





.47 


± 


0, 


,02 




J0945+2640 


B2 0942+26 





.57 


J094530 


.97+264052 


.9 


21 


.86 


± 


0, 


12 





.09 


± 


0, 


,03 




J0951+1154 


TXS 0948+121 





.37 


J095133 


.84+115459 


.6 


21, 


.81 


± 


0, 


,09 





.60 


± 


0, 


,10 




J1008+2401 


B2 1005+24 





.43 


J100832 


.62+240119 


.6 


22 


.35 


± 


0, 


,16 





,82 


± 


0, 


,20 




J1010+4159 


B3 1007+422 





42 


J101024 


.79+415933 


.3 


20, 


.69 


± 


0, 


07 





.49 


± 


0, 


,14 




J1019+4408 


B3 1016+443 





,34 


J101948 


,17+440824, 


,5 


22 


,27 


± 


0, 


,18 





,09 


± 


0, 


,05 




J1023+0424 


TXS 1021+046 





33 


J102337 


,55+042413 


,9 


21 


,62 


± 


0, 


08 





,67 


± 


0, 


,11 




,11033+3935+ 


B2 1030+39 





,41 


J103322 


,03+393551 


,0 


21 


,45 


± 


0, 


09 





,20 


± 


0, 


,02 


1.095 
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Table 1 — Continued 



Object a 


Radio Name 




SDSS Name 


SDSS r 


SDSS £phot ^spcc 






(Jy) 









7 1 no 4 i 1 1 in 

J lUo4+lllz 


T^VO 1 no 1 i 1 1 A 

LAb lUol+114 


1 

1. 


O 1 

zl 


t i no /i n c 

J 103405 


no i i i i o o 1 

,Uo+lllzol 





OO 

ZZ 


1 

Zl 


1 

± 


n 

u. 


1 *7 

1 / 


u 


o n 

39 


1 


n 

u 


1 I 

.14 






JlU4o+U5o7 


4U +1)5.45 


n 
(J 


,o7 


T 1 f\ A *y A f\ 

JlU4o4U 


on i fico^i o 

,oU+U5o71z. 


c 
,0 


o n 
ZZ 


.61 


1 

± 


L), 


1 <7 
1 / 


n 
(J 


r o 
,00 


1 


u, 


1 (■ 

lb 






J 1U45+U455 


lAb lU4o+U51 


U 




oo 


71 n/i c 1 
J 1U4551 


1 A 1 CiA K C C 1 


i 


O 1 

zl 


DZ 


1 

± 


n 
I). 


no 
Uo 


(J 


,4o 


1 

± 


u, 


1 o 

lo 






J !U4o+o45 7 


Af~\ i or no 

40 +o5.zo 


i 
1 


Uo 


71 C\AQO A 

J 1U4oo4 


O/l i O /I C70 /I 

,z4+o4o 7 z4. 


ri 

,y 


o n 
zU 


oo 
OO 


1 

± 


n 
I). 


.1)4 


(J 


on 


1 

± 


n 

u, 


nn 

,uy 


1 


KCi A 

,5y4 


T1 1 on i ooo^t 


/I i OO OT 

4U +23.27 


1. 


o o 
.OO 


711 on A o 

J HzU4o 


no i ooo^cr: 

,02+232755 


,Z 


O 1 

zl 


O 1 

.81 


1 

± 





i r 

lb 


n 

u 


58 


1 


n 




O I 

.24 


1 


o i n 

,819 


T1 1 nr i 1 nro 


t> Ty o 1 1 no i o n i 

rKb llzo+zUl 


n 

u 


1 

4o 


71 1 orCC 

J112555 


oo i 1 n C o /i 4 

.zz+195o44 


n 

.0 


O 1 

Zl 


1 O 


1 


r\ 
(.). 


1 n 

1U 


n 

u 


2 b 


1 


u. 


1 o 

lz 






J llz 7 +5 /4o 


rpYQ 1 1 O/l l (^CVn 


u 


,00 


J112/43 


, f l+5f4olo 


2 


O 1 

zl 


on 

,zy 


1 

± 


I). 


no 

Uo 


u 


53 


1 

± 


n 

u, 


no 

.Uz 






T 1 1 on i trcoo 
J 11/9+00.50 


T'vc i 1 or i crn 

IXb 1126+569 


n 




r n 

.50 


J 112904 


icr i r o o /i /i 

.15+563844 


i 

4 


O 1 

Zl 


on 

.89 


1 


0. 


o n 

20 


n 




b/ 


1 





1 n 

,10 






T1 1 A O i noocr 

J 1142+0235 


mvc i 1 on i noo 

LAb 1139+028 


0. 


o o 

38 


ti i /i onr 

J 11420b 


oo i nooroo 

.38+023533. 


.5 


o n 

20 


nn 

.90 


1 


0. 


n r 

,05 





53 


1 


0. 


n r 

,05 






J1147+4818 


A /~1 i A O O O 

4C +48.33 


0. 


r n 

.50 


J 114752 


ot i /ioiO/in 

.27+481849. 


.5 


OO 

22 


O 1 

.21 


1 


0. 


o n 

20 


0. 


1 ^ 

1 i 


1 


0. 


n r 

.05 






71 i /in i 1 ,4 n /I 

J 1148+1404 


rn "V" C 1 1 A r i 1 /I O 

IXb 1145+143 


0. 


o o 

33 


J 114825 


/io i 1 /i n /i /I n 

.43+140449 


3 


O 1 

21 


.00 


1 


0. 


11 


0. 


.40 


1 


0. 


14 






T1 OHO i 1 0(1^7 

J 1zUz+1zU7 


T'VC 1 onn i 1 O/i 

lAb 1zUU+1z4 


n 




3/ 


7 1 ono c o 

J 1zUz5z 


nn i i on *7 on 

.09+120720 





o n 

20 


b2 


1 


0. 


,04 


n 




o o 

,33 


1 

± 





1 O 

13 






J 1203+4632 1 


t~> o i onn i AHO 

133 1200+468 





1 O 

42 


t 1 ono o 1 

J 120331 


on i a r 1 o o c c 

,80+463255 


b 


OO 

22 


O O 

,23 


1 


0. 


1 O 

,13 





r n 

50 


1 


0. 


1 O 

12 






T1 OA7 i C /i r\'"7 

J 1207+5407 


4C +54.26 





.00 


T 1 OATI /I 

J 120714 


i o i c /i n 7 c /i 

,13+540754. 


5 


o n 

22 


11 


1 


0. 


18 





,77 


1 


0. 


14 






hoi c i i 7on 

J1215+1730 


4U +17.54 




n o 

03 


7 1 o i i /i 

J lz!514 


*7n i 1 to nno 

.70+173002 


,Z 


O 1 

Zl 


. 1 1 


1 


0. 


1 n 

,10 


n 




,07 


1 





n o 

03 






tiooo i cr o /i o 

J 1228+5348 


4C +54.28 





r o 

53 


J 122850 


rn i tro^om 

.60+534801 


b 


o n 

22 


r o 

,52 


1 

± 


0. 


o n 

20 





o r 

35 


1 

± 


0. 


i n 

19 






J 1238+0845 


T'VC i ooc i nnn 

IXd 1235+090 


n 

u 


41 


71 o o o 1 n 

J 123819 


or i no A rnn 

,26+084500 


4 


21 


r o 

52 


i 

± 





i n 

10 


n 




^ 5 


± 


n 




1 n 

,10 






T 1 OAfl i r-non 

J 1300+5029 


TXS 1258+507 





38 


T 1 o nn /( 1 

J 130041 


O/i i p*nonoT 

,24+502937. 





21 


,55 


± 


0. 


.07 





O" 


± 


0. 


14 


1 


,561 


J1312+1710 


TXS 1310+174 





.34 


J131235 


,21+171055. 


.7 


22 


.52 


± 


0. 


,16 





,47 


± 


0. 


,16 






J1315+0222 


TXS 1312+026 





,52 


J131516 


,93+022221 


.1 


20 


,99 


± 


0. 


,13 





,51 


± 


0. 


,17 






J1341+1032 


4C +10.36 





.69 


J134104 


,36+103207 


3 


21 


,99 


± 


0. 


,14 





.38 


± 


0. 


,11 






J1345+5846 


4C +59.20 





,42 


J134538 


,32+584654 


.6 


22 


.01 


± 


0. 


,17 





,58 


± 


0. 


.22 






J1347+1217 t 


4C +12.50 


4 


.86 


J134733 


.35+121724. 


.2 


15 


.35 


± 


0. 


00 





.18 


± 


0. 


.05 





,122 


J1348+2415 t 


4C +24.28 





.56 


J134814 


.88+241550 


.3 


22 


.31 


± 


0. 


,18 





.21 


± 


0. 


,17 


2 


,879 


J1352+3126 t 


4C +31.43 


3 


.70 


J135217 


.87+312646 


.4 


14 


.12 


± 


0. 


00 





.09 


± 


0. 


,03 





,045 


J1354+5650 


4C +57.23 





.72 


J135400 


.10+565005 


.1 


22 


.60 


± 


0. 


20 





.21 


± 


0. 


.08 






J1357+0046 


PKS 1355+01 


2 


.00 


J135753 


.71+004633 


.3 


22 


.05 


± 


0. 


,11 





.57 


± 


0. 


,18 






J1410+4850 


TXS 1408+490 





.33 


J141036 


.04+485040. 


.4 


20 


.99 


± 


0. 


,05 





.52 


± 


0. 


.09 






J1413-0312t 


TXS 1410-029 





.33 


J141314 


.87-031227. 


.3 


12 


.01 


± 


0. 


00 





.12 


± 


0. 


,06 





,006 


J1413+1509 


TXS 1411+154 





.50 


J141341 


.70+150939 


.9 


20 


.96 


± 


0. 


,09 





.43 


± 


0. 


,11 






J1414+4554 


B3 1412+461 





.41 


J141414 


.84+455448. 


.7 


20 


.16 


± 


0. 


.03 





,48 


± 


0. 


,03 





,186: 


J1415+1320t 


PKS 1413+135 


1. 


18 


J141558 


.81+132023. 


.7 


19 


.05 


± 


0. 


02 





.29 


± 


0. 


.05 





,247 
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Table 1 — Continued 



Object a 


Radio Name 


-^20cm 


SDSS Name 


SDSS r 


SDSS Zphot %i 






(Jy) 









,11421-0246+ 


4C -02.60 


0. 


.55 


J142113 


.54-024646. 


.0 


20. 


34 


± 


0. 


10 


0. 


.50 


± 


0. 


03 






J1424+1852 


4C +19.47 





.70 


J 142409 


,71+185253. 


5 


21 


.49 


± 


0. 


08 





,52 


± 


0. 


11 






J1502+3753 


B2 1500+38 





.34 


J150234 


,77+375353 


.0 


20 


.50 


± 


0. 


06 





.40 


± 


0. 


.09 






J1504+5438 


TXS 1503+548 





.38 


J150451 


,19+543839 


.7 


20 


.29 


± 


0. 


05 





53 


± 


0. 


.10 





,621 


J1504+6000 


TXS 1502+602 


1. 


.55 


J150409 


,21+600055 


.7 


20 


.00 


± 


0. 


.02 





.47 


± 


0. 


,16 


1 


,024* 


J1523+1332 


4C +13.54 





35 


J152321 


,74+133229 


2 


21 


.76 


± 


0. 


.12 





.46 


± 


0. 


.18 






J1527+3312 


B2 1525+33 





32 


J152750 


,89+331253 


.0 


22 


.04 


± 


0. 


.10 





,44 


± 


0. 


.14 






J1528-0213 


PKS 1525-020 





.47 


J152821 


,99-021319 


.0 


21 


.75 


± 


0. 


.15 





63 


± 


0. 


07 






J1548+0808 


TXS 1545+082 





.63 


J154809 


,05+080834. 


.7 


21 


.78 


± 


0. 


.14 





.65 


± 


0. 


.12 






J1551+6405 


TXS 1550+642 





.68 


J155128 


,19+640537 


,8 


22 


.10 


± 


0. 


,15 





,47 


± 


0. 


,15 






J1559+4349 


4C +43.36 





.75 


J155931 


,21+434916 


,6 


21 


.11 


± 


0. 


.07 





.42 


± 


0. 


.09 


1 


,232* 


J1604+6050 


TXS 1603+609 





.59 


J160427 


.40+605055 


.5 


21 


.44 


± 


0. 


.10 





.59 


± 


0. 


,05 






J1616+2647 


PKS 1614+26 


1. 


.41 


J161638 


.32+264701 


.3 


22 


.46 


± 


0. 


17 





.44 


± 


0. 


.03 






J1625+4134 


4C +41.32 


1. 


.72 


J162557 


.68+413440 


.8 


22 


.35 


± 


0. 


.16 





.58 


± 


0. 


.18 


2 


,550 


J1629+1342t 


4C +13.60 





.71 


J162948 


.42+134205 


.7 


22 


,34 


± 


0. 


.18 





.44 


± 


0. 


.18 






J1633+4700 


4C +47.43 





.46 


J163315 


.02+470016 


.9 


21, 


.44 


± 


0. 


.12 





.57 


± 


0. 


.19 






J1724+3852 


B2 1722+38 





.37 


J 172400 


.53+385226 


.7 


22 


.14 


± 


0. 


.12 





.11 


± 


0. 


,05 






J2203-0021 


4C -00.79 


0. 


.61 


J220358 


.31-002148 


.1 


21, 


.21 


± 


0. 


13 


0. 


.57 


± 


0. 


.10 





,729* 



a Objects noted with superscript "f" are removed from final sample because the source is : 
(1) re-evaluated morphologically to be early-type using SDSS-DR8 (J0000-1054, J1033+3935, 
J1203+4632, J1421-0246, J1629+1342); (2) nearby at z < 0.1 (J1352+3126, J1413-0312); 
(3) a well-studied source (molecular absorption system J1415+1320, alignment-effect systems 
J1120+2327 & J1348+2415, law-z CSO J1347+1217); or (5) aligned with one lobe of a double-lobe 
radio galaxy (J0834+1700). 

b Redshifts noted with superscript "*" were obtained recently by us using DIS II or TripleSpec 
at the APO 3.5m telescope (see Figures Q] & [2|). Others come from the NASA/NED database. 
Redshifts noted with ":" are uncertain. 

c Gravitational lensing system with background quasar at z = 3.2 and foreground galaxy at 
z = 0.349. 
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Table 2. NIR photometry and optical-NIR SED. 



Object J H -f^short Observation SED zg t Hubble 

Date Type a Type b 



J0003-1053 


17. 


69 


± 


0. 


.28 


16 


46 


± 


0. 


.28 


16 


19 


± 


0. 


31 


UT071213 


Q+abs 








J0134+0003 


18. 


.88 


± 





.38 


17. 


.76 


± 





.39 


16 


.91 


± 


0. 


39 


UT071213 


G 


0. 


92 


Sa 


J0249-0759 


19. 


19 


± 





.21 


18 


.30 


± 





.22 


17 


.40 


± 


0. 


.21 


UT081102 


G+Q 








J0736+2954 


19. 


62 


± 





.27 


18 


.57 


± 





.27 


17 


.60 


± 


0. 


.27 


UT081228 


G+Q 








J0747+4618 


19. 


40 


± 





.14 












17 


.95 


± 


0. 


,14 


UT091228 


Q+abs 








J0749+2129 


18. 


22 


± 





.13 












16 


,37 


± 


0. 


,13 


UT080120 


G 


0. 


,52 


SO 


J0751+2716 


19. 


.84 


± 





.18 


19 


41 


± 





.18 


18 


18 


± 


0. 


17 


UT070304 


G 


0. 


,42 


Sc 


J0759+5312 


19. 


.66 


± 





.17 


18 


.65 


± 





.17 


18 


.09 


± 


0. 


,17 


UT080218 


G+Q 








J0805+1614 


18. 


12 


± 





11 


16 


78 


± 





11 


15 


.87 


± 


0. 


,12 


UT080218 


G 


0. 


.67 


Sc 


J0807+5327 


19. 


13 


± 





.09 












17 


73 


± 


0. 


.09 


UT080120 


Q+abs 








J0824+5413 


17. 


67 


± 





.09 












15 


,81 


± 


0. 


,10 


UT080120 


G+Q 








J0839+2403 


18. 


.83 


± 





15 


18 


.60 


± 





15 


17 


63 


± 


0. 


15 


UT070304 


G+Q 








J0843+4215 


18. 


.59 


± 





.24 












16 


,51 


± 


0. 


.25 


UT080120 


G 


0. 


60 


Sb 


J0901+0304 


16. 


17 


± 





.03 


15 


18 


± 





.02 


14, 


,45 


± 


0. 


,03 


UT080120 


G 


0. 


,21 


Sc 


J0903+5012 


19. 


16 


± 





.35 


17 


.95 


± 





.35 


17 


,44 


± 


0. 


35 


UT080218 


G+Q 








J0905+4128 


19. 


46 


± 





.17 


19 


14 


± 





.17 


17, 


.48 


± 


0. 


17 


UT080218 


G+Q 








J0907+0413 


19. 


32 


± 





.33 


18 


.59 


± 





.33 


17, 


,51 


± 


0. 


33 


UT070304 


G 


1. 


10 


Sc 


J0910+2419 


18. 


04 


± 





.10 












17, 


19 


± 


0. 


,11 


UT080120 


G+Q 








J0915+1018 


20. 


.02 


± 





.23 


19 


.31 


± 





.23 


18, 


.33 


± 


0. 


23 


UT081228 


G+Q 








J0917+4725 












18 


,54 


± 





.22 


17, 


.61 


± 


0. 


22 


UT081228 


G+Q 








J0920+1753 


19. 


.00 


± 





.31 


18 


49 


± 





.31 


17, 


,53 


± 


0. 


31 


UT090116 


G 


0. 


.68 


Sc 


J0920+2714 


15. 


.91 


± 





.04 


14 


.94 


± 





.05 


14, 


.36 


± 


0. 


,06 


UT060611 


G 


0. 


23 


SO 


J0939+0304 


17. 


.80 


± 





.09 


17 


.01 


± 





.09 


16, 


.46 


± 


0. 


,10 


UT060611 


G 


0. 


,34 


SO 


J0945+2640 


19. 


33 


± 





.21 












18, 


63 


± 


0. 


,21 


UT080418 


Q 








J0951+1154 


18. 


.86 


± 





.17 


18 


.24 


± 





.17 


17, 


37 


± 


0. 


,17 


UT081228 


G 


0. 


65 


Sc 


J1008+2401 


19. 


41 


± 





.35 


18 


.76 


± 





.36 


17 


,80 


± 


0. 


36 


UT081228 


G 


0. 


62 


Sc 


J1010+4159 


17. 


.79 


± 





.12 


17 


.09 


± 





.12 


16 


.79 


± 


0. 


13 


UT060611 


G 


0. 


71 


Sc 


J1019+4408 


20. 


.05 


± 





,27 


19 


.07 


± 





,27 


17, 


.95 


± 


0. 


.27 


UT081228 


Q+abs 








J1023+0424 


18. 


.74 


± 





.15 


18 


,56 


± 





.15 


17 


63 


± 


0. 


.15 


UT060611 


G 


0. 


,71 


Sc 


J1034+1112 


20. 


.19 


± 





.29 












17 


.94 


± 


0. 


,29 


UT090116 


G+Q 








J1043+0537 


19. 


.75 


± 





.29 












17, 


.96 


± 


0. 


,29 


UT090116 


G 


0. 


42 


E 


J1045+0455 


19. 


.67 


± 





.15 


18 


.76 


± 





.15 


18, 


,53 


± 


0. 


.15 


UT081228 


Q+abs 
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Table 2 — Continued 



Object 


J 




H 






-"■short 




Observation 
Date 


SED 

Type a 




Hubble 
TvDe^ 1 


J1048+3457 


18.39 


± 


0.08 


17.74 


± 


0.13 


17.31 


± 


0.12 


UT060611 


Q+abs 






J1125+1953 


20.13 


± 


0.17 








> 19.08 


± 


0.26 


UT090116 


Q 






J1127+5743 


18.78 


± 


0.12 


17.83 


± 


0.12 


17.12 


± 


0.15 


UT060611 


G 


0.49 


Sb 


J1129+5638 


18.57 


± 


0.31 


18.02 


± 


0.31 


17.38 


± 


0.31 


UT081228 


G 


0.64 


Sa 


J1142+0235 


18.90 


± 


0.09 


17.30 


± 


0.08 


17.04 


± 


0.09 


UT060611 


G 


0.37 


Sa 


J1147+4818 


19.92 


± 


0.28 








17.23 


± 


0.28 


UT090116 


Q 






J1148+1404 


19.35 


± 


0.20 








17.72 


± 


0.20 


UT090116 


G 


0.54 


Sc 


J1202+1207 


19.36 


± 


0.08 








17.27 


± 


0.10 


UT090116 


G+Q 






J1207+5407 


18.92 


± 


0.31 


18.26 


± 


0.31 


17.58 


± 


0.31 


UT081228 


G 


0.61 


Sc 


J1215+1730 


19.78 


± 


0.16 


19.04 


± 


0.16 


18.11 


± 


0.17 


UT080120 


Q+abs 






J1228+5348 


19.31 


± 


0.36 








19.03 


± 


0.36 


UT090116 


G+Q 






J1238+0845 


17.93 


± 


0.16 


17.58 


± 


0.16 


16.88 


± 


0.17 


UT060611 


G 


0.65 


Sc 


J1300+5029 


19.27 


± 


0.13 








17.60 


± 


0.13 


UT080517 


Q 






J1312+1710 


19.55 


± 


0.37 


19.02 


± 


0.37 


18.18 


± 


0.37 


UT080218 


G 


0.63 


Sc 


J1315+0222 








16.91 


± 


0.23 


16.30 


± 


0.23 


UT060611 


G 


0.85 


Sc 


J1341+1032 


20.02 


± 


0.25 


> 19.24 


± 


0.39 


17.79 


± 


0.25 


UT081228 


Q+abs 






J1345+5846 


18.86 


± 


0.30 








18.25 


± 


0.30 


UT091228 


G+Q 






J1354+5650 


19.82 


± 


0.30 








18.14 


± 


0.30 


UT080517 


G 


0.65 


Sc 


J1357+0046 


19.77 


± 


0.19 


18.63 


± 


0.19 


18.09 


± 


0.19 


UT081228 


G+Q 






J1410+4850 








18.39 


± 


0.10 


17.08 


± 


0.10 


UT060630 


G+Q 






J1413+1509 


18.63 


± 


0.16 








17.37 


± 


0.18 


UT080120 


G 


0.22 


SO 


J1414+4554 


17.62 


± 


0.06 


16.63 


± 


0.06 


16.15 


± 


0.06 


UT070603 


G 


0.38 


SO 


J1424+1852 


19.26 


± 


0.17 


18.97 


± 


0.18 


17.77 


± 


0.19 


UT080218 


Q+abs 






J1502+3753 


18.37 


± 


0.10 


17.87 


± 


0.09 








UT060630 


Q+abs 






J1504+5438 


18.09 


± 


0.09 


17.00 


± 


0.09 


16.15 


± 


0.09 


UT080218 


G 


0.51 


Sc 


J1504+6000 


18.06 


± 


0.07 


18.46 


± 


0.05 


17.60 


± 


0.06 


UT070528 


Q 






J1523+1332 








18.51 


± 


0.21 


18.26 


± 


0.21 


UT070528 


Q+abs 






J1527+3312 














18.11 


± 


0.21 


UT080517 


G 


0.32 


Sb 


J1528-0213 








17.62 


± 


0.25 


16.75 


± 


0.25 


UT070528 


G 


0.54 


E 


J1548+0808 


18.77 


± 


0.25 


17.72 


± 


0.25 


16.83 


± 


0.25 


UT070528 


G 


0.51 


E 


J1551+6405 














18.88 


± 


0.27 


UT110518 


G+Q 






J1559+4349 


18.55 


± 


0.13 


17.70 


± 


0.12 


16.78 


± 


0.12 


UT060630 


Q+abs 
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Table 2 — Continued 



Object 


J 


H 


-^short 


Observation 


SED 


Zfit 


Hubble 










Date 


Type a 




Type b 


J1604+6050 


18.94 ± 0.18 


17.74 ± 0.18 


16.89 ± 0.18 


UT080218 


G 


0.42 


SO 


J1616+2647 




18.59 ± 0.27 


18.10 ± 0.27 


UT070528 


G 


0.72 


Sc 


J1625+4134 


19.21 ± 0.26 


17.98 ± 0.26 


16.80 ± 0.26 


UT060630 


Q 






J1633+4700 


19.10 ± 0.20 


18.49 ± 0.20 


17.73 ± 0.20 


UT080218 


Q+abs 






J1724+3852 


19.96 ± 0.09 


18.50 ± 0.09 


18.72 ± 0.09 


UT081102 


Q+abs 






J2203-0021 


18.54 ± 0.18 


17.36 ± 0.18 


16.85 ± 0.18 


UT081102 


G 


0.71 


Sc 



a The optical-NIR SED is classified as G if it drops at the blue end, Q if it keeps increasing or 
Q + abs if it flattens out at the blue end, and G + Q if there is flux excess at the red end on top of 
Q feature at the blue end. See Section [4T21 for details. 

b Sources classified as G are fit with template spectra (Figure IHJ) to estimate the Hubble type of 
the host galaxy. See section 4.2 for details. 



Table 3. New Spectroscopic Redshifts from the Apache Point 3.5m Telescope. 



Object z spec Zphot Instrument Used Lines Seen 



J0003-1053 


1.474±0.001 





,42 


TripleSpec 


[0 III], Ha, [Nil] 




T AO A T i -| / > 1 A 

J0805+1614 


0.632:±0.002 





.49 


Iripiebpec 


Ha, He I 1.08//; lines are weak & need confirmation 


J0824+5413 


0.6385±0.0003 


o 


.4(3 


DIS II 

TripleSpec 


[0 II], Call H&K, H/3, [0 III] 
Ha, [Nil], [SII], [S III]: 0.95// 




J0901+0304 


0.2872±0.0001 





.29 


DIS II 


[0 II], H/3, [0 III], [0 I], Ha, [N II], [S II] 
also detected in H I @z=(0.2869, 0.2886); see 


Paper 3 


J0920+2714 


0.2064±0.0002 





.17 


DIS II 


Ca II H&K, G-band, Ha, [N II], [S II] 
also detected in H I @z=0.2067 (Paper 3) 




J1504+6000 


1.024±0.001 





.67 


DIS II 


C III], C II]2326A, [Ne V], [0 II], [Ne III] 3869 
slightly revises Burbidge (1970) redshift of z- 


, 3967A; 
=1.022 


J1559+4349 


1.232±0.001 





.42 


TripleSpec 


[0 III], Ha, [N II], [S III] 0.95//, He I 1.08//: 




J2203-0021 


0.729±0.001 





.57 


TripleSpec 


Ha, [Nil], [S III] 0.95// 
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Table 4. Sources with Discrepant Redshift Estimates. 



Object 




Hubble 


TC c 
-"■ short 


2K 


Az 






Type b 








J1008+2401 


0.62 


Sc 


17.80 


1.25 


0.63 


J1043+0537 


0.42 


E 


17.96 


1.35 


0.93 


J1142+0235 


0.37 


Sa 


17.04 


0.85 


0.48 


J1148+1404 


0.54 


Sc 


17.72 


1.20 


0.66 


J1312+1710 


0.63 


Sc 


18.18 


1.52 


0.89 


J1354+5650 


0.65 


Sc 


18.14 


1.49 


0.84 


J1413+1509 


0.22 


SO 


17.37 


1.00 


0.78 


J1527+3312 


0.32 


Sb 


18.11 


1.46 


1.14 


J1616+2647 


0.72 


Sc 


18.10 


1.46 


0.74 



a ' 6 ' c Same as in Table 2. 
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Fig. 1. — The APO/DIS optical spectra with prominent emission/absorption lines marked. 
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Fig. 2. — The APO/TripleSPEC near-IR spectra with prominent emission lines marked. The 
shaded regions are regions of significant atmospheric absorption. 
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SDSS composite /-band A" s -band 




Fig. 3. — Representative optical and NIR images from our sample all shown at the same scale and 
orientation as indicated in left-hand column. Shown are the SDSS composite images (left; blue: 
g-band, green: r-band, red: z-band), and the APO 3.5m J-band (middle) and K s band (right) 
images of (from top to bottom) J0805+1614, J0901+0304, J0920+2714, and J0939+0304. 
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Fig. 4.— SED classified as Q {top, J1125+1953) and Q + abs (bottom, J1724+3852). 
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Fig. 5.— SED classified as G {top, J0951+1154; bottom, J0920+2714). The spectra are template 
galaxies of Sc type at z = 0.65 (top) and SO type at z = 0.23 (bottom). 
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Fig. 6.— SED classified as G + Q (J0917+4725). 
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Fig. 7. — Top (a): Spectroscopic redshifts compared with photometric redshifts from DR8 for 
those 16 sources with spectroscopic redshifts. The line shows z spcc = z p h ot , illustrating that the 
fitting algorithm does not work for high-z {z > 1) quasars. Even for those sources not classified as 
quasars, the z p hot is systematically lower than z spec by Az = 0.05-0.1 at z ~ 0.5-1. Bottom (b): 
Photometric redshifts from DR8 (asterisks with error bars) and from our fitting routine (diamonds) 
using SDSS+NIR photometry for sources with spectroscopic redshifts and classified as G. See 
Section 4.3 for discussion. 
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Fig. 8. — Template spectra for typical galaxies (Mannucci et al. 2001) and unreddeded QSOs 
(Vanden et al. 2001, Glikman et al. 2006). 




Fig. 9. — The observed K s — z diagram for objects whose SED is classified as G. The red dots 
represent objects with known spectroscopic redshifts. For the remainder, photometric redshifts are 
used (black asterisks). The black line is the best fit K — z relation from Willott et al. (2003) with 
the dotted lines indicating the la dispersion and the dashed lines indication the 2a dispersion. 
Table 4 lists 9 sources which are fainter in K than the 2a dispersion in the (K — z) relation. See 
Section 4.3 for discussion. 



